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EXECUTIVE  SUMMARY 


The  ICDA  (Initial  Capability  Decision  Aid)  is  a  heat  stress  prediction  model 
recently  developed  for  monitoring  the  physiological  status  of  Soldiers.  The  ICDA  is  an 
abridged  version  of  the  USARIEM  human  thermal  physiological  simulation  computer 
model  SCENARIO-J  (Kraning  and  Gonzalez.,  1997)  and  the  Gagge  model  (Gagge  et 
al.,  1986).  ICDA  utilizes  real  time  input  of  metabolic  activity  (M)  derived  from  heart  rate 
(HR)  and  air  temperature.  By  using  HR  to  estimate  M,  Soldiers  required  to  perform 
multiple  tasks  over  extended  time  periods  can  be  monitored  with  non-invasive  sensors 
which  are  both  convenient  to  place  and  comfortable  to  wear.  This  study  compares  ICDA 
predictions  of  core  temperature  (Tcr),  skin  temperature  (Tsk),  and  sweat  rates  (SR)  with 
laboratory  measured  values. 

The  model  validation  was  conducted  using  data  from  three  laboratory  studies 
with  varied  environments  (i.e.,  27°C/75%RH;  49°C/18%;  35°C/45%),  clothing 
configurations  and  heat  acclimation  status.  Subjects  walked  on  a  treadmill  at  -270-628 
W  for  up  to  180  min,  depending  upon  study  conditions.  Tcr  rise  was  restricted  to 
<39.5°C.  Individual  anthropometries,  physiological  and  environmental  time  series  data 
were  collected  in  each  study.  Root  Mean  Square  Deviation  (RMSD)  was  utilized  to 
examine  the  difference  between  laboratory  measurements  and  ICDA  predictions  of  Tcr 
and  Tsk  across  time.  In  addition,  predicted  SR  was  compared  to  the  measured  SR 
where  available. 

RSMD  of  Tcr  ranged  between  ±0.1 1  and  ±0.50°C  for  these  studies.  The  model 
predicted  measured  values  of  Tcr  within  one  standard  deviation,  when  environmental 
conditions  were  less  stressful.  Errors  of  estimates  were  larger  with  small  sample  size 
(N<5).  The  ICDA  predictions  of  Tsk  were  less  accurate  than  those  of  Tcr.  Generally, 
model  predictions  of  Tsk  were  more  difficult  than  Tcr  because  measured  Tsk  varied  by 
skin  location,  clothing,  environmental,  and  operational  conditions.  ICDA  predictions  of 
SR  were  more  accurate  when  environmental  conditions  were,  in  terms  of  heat  stress, 
less  extreme.  However,  for  encapsulated  subjects,  working  in  a  35°C/45%RH 
environment,  the  results  for  predictions  of  SR  were  mixed. 

Predictions  of  thermal  status  of  Soldiers  during  heat  stress,  particularly  for  Tcr, 
using  this  simplified  model,  are  promising.  To  improve  predictions  of  Tsk,  the  existing 
data  bases  should  be  evaluated  to  determine  which  Tsk  measurement  site  best 
represents  the  combined  Impact  of  individual  variability,  clothing,  environment,  and 
operations  on  Tsk,  and  is  therefore,  the  best  site  to  use  for  model  development  and 
validation  of  Tsk  predictions.  Similarly,  for  SR  predictions,  the  trends  of  individual  SR 
patterns  need  to  be  carefully  evaluated  to  further  refine  this  model.  Lastly,  it  is 
recommended  that  ICDA  be  applied  to  real  time  situations  in  both  laboratory  and  field 
training  exercises  for  future  Tcr  sensitivity  analysis,  and  further  model  development. 
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INTRODUCTION 


The  ICDA  (Initial  Capability  Decision  Aid)  model  is  a  heat  stress  prediction  model 
recently  developed  for  monitoring  the  physiological  status  of  Warfighters.  It  is  a  basic 
and  abridged  model  derived  from  elements  of  the  USARIEM  human  thermal 
physiological  simulation  computer  model,  SCENARIO-J  (Kraning  and  Gonzalez,  1997), 
and  the  Gagge  Model  (Gagge  et  al.,  1986).  The  model  can  be  loaded  with 
anthropological  characteristics  of  an  individual  (i.e.,  height,  weight,  and  clothing),  group 
means  or  use  default  population  values.  Real-time  inputs  of  measured  HR  and  local 
weather  (e.g.,  temperatures,  wind,  relative  humidity,  estimates  of  radiant  load)  are  used 
by  the  computer  model  to  make  real-time  predictions  and  estimates  of  subjects’ 
physiological  status  or  other  physiological  parameters. 

A  main  purpose  for  constructing  this  model  was  to  predict  time  estimates  of  the 
internal  body  temperature  (Tcr),  sweat  rates  (SR),  and  hydration  status  of  Soldiers  and 
Warfighters  in  a  battlefield  situation.  For  instance,  Tcr  is  the  traditional  and  common 
physiological  parameter  for  heat  strain  assessment  as  it  reliably  indicates  impeding 
injury  (Amos  et  al.,  2000;  Pandolf  and  Goldman  1978).  However,  measuring  Tcr  is 
invasive  and  may  be  impractical  for  real  time  monitoring  of  Soldiers  engaged  in  long 
hours  of  various  multiple  and  unpredictable  tasks  in  widely  disbursed  hot  field 
environment.  Although  obtaining  “true”  values  of  many  parameters  from  environmental, 
physiological  and  operational  conditions  are  desirable  and  increases  the  model’s 
accuracy  to  assess  Soldiers’  status,  alternatives  for  the  worst  case  scenario  (e.g.,  losing 
signals  from  sensor  devices,  wrong  calibration,  conflicts  between  sensors,  etc)  still  need 
to  be  considered.  For  these  reasons,  ICDA  was  developed  to  predict  physiological 
responses  to  battlefield  situations  from  a  minimum  number  of  non-invasive  inputs. 

The  thermal  physiology  of  the  ICDA  model  is  represented  schematically  in  Figure 
1 ,  where  the  human  is  modeled  as  two  physiological  compartments  (core  and  skin) 
surrounded  by  a  passive  clothing  compartment.  Within  a  compartment  the  properties 
are  uniform:  i.e.,  everywhere  in  core  at  time  t,  Tcr  has  the  same  value.  All  metabolic 
heat  production  (M)  occurs  in  the  core.  Some  heat  is  lost  directly  from  the  core  to  the 
environment  by  respiration;  all  of  the  other  heat  from  the  core  is  transferred  to  the  skin 
by  conduction  or  convection  with  skin  blood  flow.  The  primary  method  used  to  regulate 
Tcr  is  achieved  by  controlling  blood  flow. 
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Figure  1 .  Schematic  of  core,  skin  and  clothing  compartments 


Evaporation  from  skin 


Respiratory  heat  loss 


Skin  blood  flow 


Convection  to  air 


Conduction  through  clothing 


''A 

Radiation  to  surrounding  surfaces 


=  core 


I  I  =  skin  compartment 
SSi  =  clothing  property 


The  basic  ICDA  operational  structure  is  summarized  in  Figure  2.  The  outline  identifies 
the  simulation  routines  used  to  predict  an  individual’s  Tcr,  sweat  rate,  accumulated  water 
loss  (WL),  skin  temperature  (Tsk),  metabolism,  and  a  thermoregulatory  strain  indicator 
quantified  by  the  physiological  strain  index  (PSI)  (Moran  et  al.,  1988). 
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Figure  2.  Basic  operational  structure  of  ICDA 


Recent  reports  from  wars  in  Iraq  and  Afghanistan  indicate  that  faster  access  to 
medical/surgical  care  and  identifying  evacuation  situations  would  save  more  Soldiers’ 
lives  (Bilski  et  al. ,  2003;  Gawande,  2004).  The  current  strategies  to  shorten  the  time  lag 
between  injury  and  treatment  is  to  position  a  small  surgical  care  unit  near  the  battlefield 
(Bilski  et  al.,  2003).  The  model  can  contribute  this  effort  by  forecasting  both  real-time 
and  future  probabilities  of  Soldiers’  health  status  by  notifying  unit  medical  personnel  and 
command  elements  on  the  battlefields.  The  ICDA  will  also  provide  a  longitudinal 
projection  of  SR,  which  can  be  used  to  predict  Soldiers’  dehydration  status  and  water 
requirements.  SR  predictions  are  important  for  both  logistics  and  for  Soldiers  to  assess 
unit  and  individual  water  requirements  (Montain  et  al.,  2005).  Thus,  logistics  personnel 
can  plan  for  water  delivery  to  meet  the  projected  demand  and  the  associated  cost,  and 
Soldiers  can  sustain  their  physical  and  cognitive  performance  by  preventing  from 
dehydration. 

This  study  described  the  basic  concept  of  the  ICDA  model,  and  presents  an 
evaluation  of  model  performance  with  a  focus  on  Tcr,  Tsk  and  SR  predictions  using 
available  data  from  military  heat  studies.  Comparisons  of  measured  physiological  data 
with  the  model  predictions  provide  the  analytical  basis  needed  to  characterize  model 
performance  and,  if  necessary,  identify  further  improvements  of  the  model. 
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METHODS 


ICDA  MODEL 

The  ICDA  model  as  described  in  the  Introduction  represents  the  human  by  two  active 
physiological  compartments  (core  and  skin)  surrounded  by  a  passive  clothing 
compartment.  The  model  was  designed  for  military  personnel,  who  routinely  are  heat 
acclimatized  and  have  normal  or  low  levels  of  body  fat  relative  to  the  general 
population. 

Metabolism  (M):  The  source  of  heat  production  in  the  core  compartment  is  metabolic 
activity.  Thus,  M  is  an  important  parameter  for  estimating  Tcr  and  Tsk.  Real  time 
estimates  of  M  for  input  into  the  ICDA  model  are  derived  from  the  Soldiers’  measured 
HR  and  environmental  temperatures  (Berglund,  1977)  using  the  following  equation: 

M  =  [0.68+4.69-  (HRratio-1)  -  0.052-  (HRratio-1)  •  (Ta-20)]  -58.1-Adu  [W]  (1) 

where  HRratio  =  observed  HR  given  at  the  time  divided  by  the  resting  HR  of  the 
individual,  and  Ta  =  ambient  temperature  in  °C.  Adu  is  the  body  surface  area. 
The  basis  for  Equation  1  is  that  M  as  HR  increases  for  a  given  task  at  an 
elevated  Ta.  The  heart  pumps  blood  to  transport  oxygen  to  the  cells  to  support 
the  metabolic  effort  and  to  exchange  heat  at  the  skin  in  order  to  cool  the  core. 
This  increased  HR  to  the  skin  is  a  non-metabolic  activity.  The  metabolic 
adjustment  for  Ta  in  Equation  1  corrects  HR  for  the  non-metabolic  cooling 
function  of  skin  blood  flow. 

This  equation  is  applicable  for  the  conditions  including  20°C  <  Ta  <  40°C,  1 .2  < 
HRratio  <  2.1,  wind  speed  =  1 .25  ms'1,  dew  point  temperature  <  20°C,  and  an  intrinsic 
clothing  insulation  of  about  0.6-0.7clo  (i.e.  that  of  a  BDU).  Within  these  limits,  a  good 
relationship  exists  between  HR  and  02  uptakes  for  different  work  rates  in  laboratory 
studies  (Berglund,  1977).  HR  is  non-invasive  and  convenient  to  measure  compared  to 
Tcr,  particularly  for  long-duration  duty  assignments. 

Heat  balance:  A  heat  balance  analysis  of  the  core  compartment  yields, 

M/Adu=qres+qk+qskbf+(Wc/Adu)  -cbt-  (dTcr/dt)  [W/m2]  (2) 

where 

Adu  =  the  body  surface  area 
qres  =respiratory  heat  loss 

Specific  heat  of  body  tissue  (cbt)  is  constant  as  .97W-h/°C-kg,  and 
Wc  =  the  weight  (kg)  of  core. 

Passive  heat  conduction  (qk)  from  core  to  the  skin  =k®(TcrTSk) 

Conductance  (k)  of  tissue  between  core  and  skin  =  5.28W/(°C-m2). 
the  heat  transported  by  blood  flow  (qskbf)  to  the  skin  =Skbf-cpb-(TCr-Tsk), 

Where  specific  heat  of  blood  (cpb)  =  1 .163  W-h/(L-°C). 
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Skin  blood  flow  (Skbf)  =  [Skbfn+Cdil*  (Tcr-Tcset)]/(1+Cstr»  (TSkserTSk)) 
where  The  vasolidation  coefficient  (Cdil)  =  50L/(h*m2»°C) 

The  neutral  Skbfn  =  6.3  L/(h*m2), 
vasoconstriction  coefficient,  Cstr,  =  0.5  °C'1. 

Tcset  =  36.8°C,  and  Tskset  =  33.7°C. 

Skbf  is  modeled  proportional  to  changes  in  Tcr  and  Tsk  from  set  point  temperatures 
(Tcrset, Tskset)-  The  maximum  and  minimum  Skbf  limits  were  set  to  90  and  2  L/min, 
respectively. 

The  rate  of  Tcr  change  (Tcr/dt)  found  by  rearranging  Equation  2: 

Tcr/dt=(M/Adu-qres-qk-qskbf)/(Wc*cbt/Adu),  (3) 

can  be  step-wise  integrated  to  find  the  next  Tcr  (Tcr2)  after  time  step  (At): 

TCr2  =  Tcn+fdTcr/dt]!  At  (4) 

In  a  similar  fashion,  the  energy  balances  of  the  skin  compartment  can  result  in: 

qskbf+qk=qdry+qevap+Wsk/Adu«cbt»(dTskc/dt)  (5) 

where  qskbf  and  qk  represent  conduction  and  blood  flows  from  the  core  to  skin. 
Dry  heat  flow  from  the  skin  to  the  clothing  compartment  is  qdry,  and  qevap  is  the 
evaporative  heat  losses  from  the  skin.  Wsk  represents  the  weight  (kg)  of  skin. 

As  with  the  core  compartment,  the  rate  of  change  of  skin  temperature  in  the 
compartment  determined  from  the  energy  balance  can  be  stepwise  integrated  to  find 
the  compartment’s  skin  temperature  at  time  t+At. 

Heat  transfer:  The  model  of  heat  exchange  between  the  human  skin  surface,  clothing 
and  environment  is  defined  by  functions  for  effective  air  movement,  resistance  to  heat 
flow  by  radiation  and  convection,  and  the  resistance  to  evaporative  heat  transfer 
(Kraning  and  Gonzalez,  1991).  This  dry  heat  ioss  (Dry)  was  determined  by: 

Dry  =  (Tsk-T0)/(Rcit)  [W/m2]  (6) 

where  Rcit  is  the  total  dry  thermal  resistance  between  skin  and  the  environment 
and  T0  is  the  operative  temperature.  Rdt  values  used  in  the  ICDA  validation 
process  are  from  thermal  manikin  measurements  for  the  specific  clothing 
ensembles  used  in  this  model.  T0  is  defined  as: 

T0  =  (hc«Ta  +  hr«Tr)/(hc+hr)  [°C]  (7) 
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where  Ta  and  Tr  in  Equation  7  represent  air  and  radiant  temperatures, 
respectively.  The  hr  and  he  terms  are  the  coefficients  of  radiant  and  convective 
heat  transfer,  defined  by  the  following  two  equations: 


hr  =  4cr«(Ar/AD)« 
he  =  8.6*V53 


Tel  +  To 


+  273 

2i 


[W/(°Om  )] 


[W/(°C*m2)] 


(8) 

(9) 


where  a  is  the  Stefan-Boltzman  constant  known  as  5.67X1  O'8  [W/(m2»K4)],  and 
Ar/AD  is  the  of  body  surface  fraction  exposed  to  radiation  equal  to  about  0.725  for 
a  standing  person.  For  low  wind  (V  in  m/s)  conditions,  he  is  estimated  from 
metabolic  activity  (see  Appendix  A  for  details). 


The  total  respiratory  heat  loss  is  the  sum  of  the  evaporative  respiratory  loss  (Eres) 
(Fanger,  1972)  and  convective  respiratory  loss  (Cres).  They  are  related  to  exercise 
intensity  at  a  given  environment  (Kraning  and  Gonzalez,  1991): 

Eres  +  cres  =  0.0023«M«  (44-Pa)  +  0.001 4*M«  (34-Ta)  [W/m2]  (10) 

where  M  is  total  metabolic  activity  in  watts  (W) ;  Pa  is  ambient  vapor  pressure; 
and  Ta  is  air  temperature. 


Quantitatively,  relative  to  respiratory  heat  loss,  heat  loss  from  the  skin  surface  is 
much  more  significant.  The  maximum  rate  of  evaporative  heat  loss  (Emax)  from  the  skin 
surface  is  calculated  based  on  the  vapor  pressure  difference  between  skin  surface  and 
air: 


Emax  ~  (Psk"Pa)/F?pclt  [W/m  ]  (11) 

where  Psk  =  saturated  vapor  pressure  (Torr)  of  water  at  skin  temperature;  Pa  = 
ambient  vapor  pressure;  and  Rpcit  =  total  vapor  resistance  of  the  clothing  from 
skin  to  ambient.  RpC|t  values  used  in  this  validation  study  report  were  measured 
by  a  sweating  thermal  manikin  for  the  specific  clothing  ensembles  used  in  this 
model. 


Sweat  rates  (SR):  The  sweat  rate  (SR)  for  evaporative  heat  loss  is  modeled  as  being 
proportional  to  changes  in  mean  body  ternperature(Tmb)  and  skin  temperature  from  set 
point  values: 


SR=1 70  (Tmb-Tmbset)e'((Tsk‘Tskset)/10'7)  g/(h«m2)  [g/min]  (12) 

where  the  set  point  of  Tmb,  Tmbset  =  36.49°C  and  Tskset  =  33.7°C 

Evaporative  heat  loss  (Esk)  is  calculated  from  skin  wettedness  (w),  the  fraction  of  the 
skin  covered  with  sweat,  and  the  maximum  evaporation  rate  (Emax)  of  100%  wet  skin: 
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Esk=w»Emax,  where  w=SR«hfg/Emax  and  hfg  is  latent  heat  of  evaporation  (=.68  W»h/g).  If 
SR*hfg>Emax,  w=1  and  excess  sweat  drips  off  the  skin.  For  the  validation  modeling,  the 
maximum  sweat  rate  was  limited  to  667  g/(h»m2). 

PSI:  In  addition  to  Tcr,  TSk,  and  SR,  the  ICDA  also  calculates  the  Physiological  Strain 
Index  (PSI).  PSI  combines  HR  and  Tcr  values  at  any  given  time  to  provide  a  simple, 
concise  and  ordinal  physiological  strain  scale  between  0  and  10.  PSI  has  been  useful 
to  evaluate  heat  strain  in  relation  to  various  combinations  of  environmental  (e.g.,  heat, 
humidity),  operational  (e.g.,  clothing,  work  rates)  and  biological  (i.e.,  gender)  conditions 
(Moran  et  al.,  1988).  The  PSI  was  calculated  as  follows  (13): 

PSI  =  5(Tcrt-  T cr0)*(39 .5-T cro)"1  +5*  (HR,-HR0)  •  (1 80-HRq)'1  (1 3) 

where  Tcro  and  HR0  are  the  initial  Tcr  and  HR  measurements  at  rest,  and  Tcrt  and 

HRt  are  simultaneous  measurements  at  a  given  time  (f). 

The  details  of  the  physiological,  thermodynamic  and  heat  transfer  functions  of  the 
ICDA  are  described  further  in  APPENDIX  A. 


VALIDATION  DATA 

The  model  validation  was  conducted  using  three  laboratory  heat  studies.  The 
investigators  adhered  to  the  policies  for  the  protection  of  human  volunteers  as 
prescribed  in  Army  Regulation  70-25,  and  the  research  was  conducted  in  adherence 
with  the  provisions  of  45  CFR  Part  46.  Test  volunteers  provided  their  informed  consent 
prior  to  participation  in  the  study  and  were  free  to  withdraw  from  the  study  at  any  time. 

Heat  Studyl  (HS1):  Nine  volunteers  (8  males,  1  female;  age:  23  ±4  [SD]  yr; 
height  174.2  ±  5.2  cm;  weight:  73.4  ±  6.5  kg;  Body  Mass  Index  (BMI):  24.2  ±  1.6), 
wearing  hot  weather  battle  dress  uniform  (HWBDU)  did  intermittent  exercise  in  a  warm- 
humid  environment  (27°C,75%RH)  for  170  minutes  (Santee  et  al.,  2005) ,  The  exercise 
routine  consisted  of  walking  on  a  level  treadmill  at  1 .34  ms'1  for  30  minutes  followed  by 
1 0  minutes  of  rest.  These  individuals  were  not  heat  acclimated  prior  to  the  study.  The 
study  conditions  were  relatively  benign  and  no  subjects  were  withdrawn  due  to  pre-set 
physiological  safety  limits  (Santee  et  al.,  2005).  Tcr  was  measured  rectally  and  the  skin 
surface  temperature  was  averaged  from  3  regions  (back,  arm,  thigh)  measured  by  heat 
flow  discs  (Santee  et  al.,  2005).  The  average  Tsk  was  calculated  utilizing  Burton’s 
weighting  method  (0.5»Tback+0.36«Tthigh+0.14«Tarm)  (Santee  et  al.,  2005). 

Heat  Study2  (HS2):  Seven  male  Soldiers,  wearing  Army  Physical  Fitness 
Uniform  (APFU)  (i.e.,  T-shirt,  shorts,  socks,  shoes),  continuously  walked  on  a  treadmill 
with  a  4%  grade  at  1 .56  m-s'1.  Metabolisms  ranged  from  420  -  620  watts  in  a  hot  dry 
(49°C,18%RH)  condition  (Montain  et  al.,  2005).  The  maximum  exercise  duration  was 
100  minutes  and  subjects  discontinued  a  test  at  a  point  of  their  voluntary  exhaustion  or 
Tcr  >  39.0°C.  Subjects  (age:  22  ±  5  yr,  height  176  ±  4  cm;  weight:  71.1  ±  9.3  kg;  BMI: 
22.9  ±  2.4)  were  un-acclimated  on  the  1st  day  of  the  study.  The  experiment  was 
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repeated  daily  for  10  days  (Montain  et  al. ,  2005).  Only  five  volunteers  participated  on 
the  10th  day.  The  data,  collected  on  the  1st  day  (non-acclimated  stage)  and  10th  day 
(acclimated  stage),  were  utilized  in  this  study.  An  ingestible  telemetry  temperature  pill 
was  utilized  to  measure  Tcr.  Skin  temperatures  were  measured  by  thermocouples  from 
four  regions  (i.e.,  chest,  forearm,  thigh,  and  calf)  to  calculate  a  mean  weighted  TSk 
based  on  Ramanathan  (1964)  method  described  as  0.3*  (TChest+Tarm)  +  0.2*  (Tthigh+Tieg) 
(Montain  et  al.,  2005). 

Heat  Study3  (HS3):  Eight  heat  acclimated  men  (age:  23  ±  6yr;  height:  176  ±  6 
cm;  weight:  76.0  ±15.4  kg;  BMI:  24.4  ±  4.3),  wearing  chemical  protective  garments, 
walked  on  a  treadmill  with  a  4-9  %  grade  at  1 .56-1  .eSrms'1.  Metabolic  effort  averaged 
450W  and  oxygen  uptake  was  -55%  of  Vo2max.  Participants  walked  until  they  voluntary 
withdrew  from  the  study  or  Tcr  >  39.5°C  (Latzka  et  al.,  1998).  Subjects  were  able  to 
continue  exercising  for  <  40  minutes  at  35°C/45%RH.  Their  physical  characteristics 
(mean  ±  SD)  including  age  and  anthropometry  were:  23  ±  6  yr,  height:  174.2  ±  5.2  cm; 
weight:  73.4  ±  6.5  kg  (Latzka  et  al.,  1998).  Tcr  was  measured  rectally.  Tsi<  was 
measured  by  thermocouples  at  four  sites  (i.e.,  forearm,  chest,  thigh,  calf)  and  mean  TSk 
was  calculated  based  on  Ramanathan  (1964)  weighting  method  (Latzka  et  al.,  2005). 

Table  1  shows  the  summary  of  these  three  studies.  Depending  upon  the  study 
schemes,  HR,  Tcr,  and  TSk  were  collected  at  different  time  intervals.  Volunteers  in  all 
studies  exercised  without  water  replacement.  Sweat  rates  were  measured  in  HS1  and 
HS3,  calculated  by  the  difference  in  body  weight  and  associated  weights  (i.e.,  clothing) 
between  before  and  after  work.  If  resting  HR  was  not  available,  a  default  value  of  70 
bpm  was  used  as  the  model  input;  otherwise  actual  HR  of  subjects  who  were  resting  or 
sitting  at  the  beginning  of  the  studies  or  the  lowest  HR  during  the  studies  was  utilized  as 
resting  HR.  In  addition,  if  measured  HR  values  during  more  active  phases  were 
missing,  values  interpolated  over  the  time  interval  between  existing  data  points  were 
substituted.  The  HR  summary  in  each  study  is  displayed  in  APPENDIX  B, 
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Table  1 .  The  data  summary  used  in  this  study 


Study# 

HS1 

HS2 

HS3 

Ta  (°C) 

27 

49 

35 

RH  (%) 

75 

18 

45 

Acclimation  Status 

No 

No->Yes 

Yes 

Clothing 

HWBDU 

T-shirt  &  Shorts 

Protective  clothing 

Activity  (W) 

277-350 

412-628 

343-552 

Duration  (min) 

-170 

<100 

<40 

n 

9 

7 

8 

Source 

Santee  et  al.,  (2005) 

Montain  et  al.,  (2005) 

Latzka  et  al.,  (1998) 

STATISTICAL  ANALYSIS 


The  predicted  and  observed  Tcr  and  TSk  were  compared  using  Root  Mean  Square 
Deviation  (RMSD)  in  each  individual.  The  RMSD  was  used  to  quantify  the  average 
difference  between  predicted  and  observed  measurements  across  time  (4). 


The  RMSD  was  calculated  as  follows: 


RMSD  = 


(14) 


where  dj=  difference  between  observed  and  predicted  at  1-min  intervals;  and  n  = 
the  number  of  compared  points. 

In  addition,  Evaporative  Sweat  Rates  (EvapSR)  predicted  by  the  ICDA  model  were 
compared  with  measured  EvapSR  in  available  studies,  using  a  paired  t-tests. 
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RESULTS 


CORE  TEMPERATURE 

Figure  3  (a-d)  summarizes  the  comparisons  of  mean  measured  Tcr  to  the 
corresponding  mean  ICDA  predictions.  RMSD  ranged  between  0.1 1  and  0.50  °C  for 
these  studies.  The  model  predictions  were  within  one  standard  deviation  of  the  mean 
measured  values,  when  Tcr  of  the  volunteers,  wearing  BDU  and  walking  at  27°C  /75% 
environmental  conditions,  were  less  than  38.0°C  (e.g.,  Figure  3a).  For  the  HS2  data, 
prediction  errors  were  greater  in  non-acclimation  subjects  than  in  acclimated  subjects  at 
49°C/18%  (e.g.,  Figure  3b  and  3c).  Although  Tcr  in  unacclimated  individuals  in  HS2 
quickly  reached  very  high  levels  (>  39.0°C)  (Figure  3b),  after  the  10-day  heat 
acclimation  process,  the  same  individuals  were  able  to  maintain  their  Tcr  at  a  lower  level 
(~38.5°C)  (Figure  3c).  Initially,  for  heat  acclimated  Soldiers  at  49°C/18%RH,  the  model 
predicted  a  rate  of  increase  in  Tcr  greater  than  that  measured;  however,  the  model  was 
more  accurate  toward  the  end  of  exercise  (Figure  3c). 

Errors  of  estimates  were  also  larger  when  sample  sizes,  due  to  the  voluntarily 
terminations,  were  smaller  (n<5)  (e.g.,  Figure  3b,  3d).  However,  ICDA  predictions  were 
reasonable  for  HS3  in  which  subjects,  wearing  protective  garments  and  walking  in  a  hot 
humid  condition,  increased  their  Tcr  <  39.5°C. 
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Figure  3  The  summary  comparisons  between  mean  measured  and  predicted  core 
temperature  (Tcr)  from  different  heat  studies. 
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Figure  3  cont. 

c.  Heat  Study  2  (daylO  -  heat  acclimated  state;  49°C,  18%;  army  physical  fitness 
uniform;  n=5) 
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d.  Heat  Study  3  (35°C,  45%;  protective  garment;  n<8) 
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SKIN  TEMPERATURE 


The  results  of  ICDA  predictions  for  Tsk  varied  between  the  different  studies.  Tsk 
measurements  on  the  1st  day  in  HS2  were  not  available.  Overall,  in  comparison  to  the 
Tcr  predictions,  the  model  predictions  of  Tsk  were  less  accurate.  Figure  4  (a-c)  shows 
the  comparisons  between  measured  Tsk  and  the  model  predictions  for  the  three  studies. 
In  all  studies,  the  model  tended  to  over-predict  the  measured  mean  weighted  Tsk  (Figure 
4a-c).  In  comparison  to  the  other  studies,  predictions  for  HS3  (Figure  4c)  showed  better 
agreement  with  measured  mean  weighted  Tsk(RSMD  =  0.64°C).  During  HS3, 
volunteers  wore  protective  clothing  which  limited  evaporative  heat  transfer  through 
clothing  to  the  environment.  In  contrast,  for  HS1  and  HS2  (Figure  4a,  4b),  the  model 
tended  to  over-predict  measured  mean  weighted  Tsk  by  about  1  °C  or  more. 

Figure  4.  Comparisons  between  mean  measured  and  predicted  skin  temperature  (Tsk) 
for  three  heat  studies. 

a.  Heat  Study  1  (27°C,  75%;  hot  weather  battle  dress  uniform;  n=9) 
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Figure  4  cont. 

b.  Heat  Study  2  (daylO  -  heat  acclimated  state;  49°C,18%;  army  physical  fitness 
uniform;  n=5) 
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c.  Heat  Study  3  (35°C,  45%;  protective  garment;  n<8) 
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Accurate  predictions  for  Tsk  were  problematic  because  measured  Tsk  varied  by 
skin  locations  and  between  individuals.  Figure  5a  and  Figure  5b  are  examples  of 
regional  Tsk  collected  from  two  volunteers  in  HS1 .  The  first  individual’s  thigh 
temperature  (Tthigh)  (Figure  5a)  agreed  closely  with  predicted  Tsk.  However,  the  back 
(Tback)  and  arm  (Tarm)  temperatures  were  cooler  after  30  minutes  and  the  deviation 
between  observed  and  predicted  values  for  the  arm  increased  with  time.  The  second 
subject’s  T ,hjgh  (Figure  5b)  also  agreed  well  with  the  predicted  Tsk.  The  TbaCk  was  similar 
to  Tarm  until  the  Tarm  sensor  became  detached  during  the  final  140-160  minute. 

Figure  5a.  The  distributions  of  measured  skin  temperatures  by  different  regions  and 
ICDA  predictions  in  one  individual  (SN7)  from  HS1:  Example  1. 
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Figure  5  cont. 

Figure  5b.  The  distributions  of  measured  skin  temperatures  by  different  regions  and 
ICDA  predictions  in  one  individual  (SN4)  from  HS1 :  Example  2. 

- ICDA-Tsk4  *■  -  Tback4  — —  Tarm4  Tthigh4  ▲  wtedTsk4 


36 


0  20  40  60  80  100  120  140  160  180 

minutes 

As  opposed  to  the  various  Tsi<  responses  in  HS1  (Figure  5a,  b),  the  responses  of 
TSk  by  region  were  relatively  uniform  in  HS3  when  subjects,  wore  protective  garments. 

Figure  6a  and  6b  are  examples  of  measured  Tsk  by  region  in  HS3.  Both 
individuals  showed  lower  Tcaif  and  higher  Tchest  than  any  other  Tsk  region.  The  maximum 
difference  between  Tcaif  and  TChestcan  be  4°C;  however,  the  regional  Tsk  responses  to 
heat  stress  during  the  experiment  were  very  similar. 
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Figure  6a.  The  distributions  of  measured  skin  temperatures  by  different  regions  and 
ICDA  predictions  in  one  individual  (SN2)  from  HS3:  Example  1 . 
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Figure  6b.  The  distributions  of  measured  skin  temperatures  by  different  regions  and 
ICDA  predictions  in  one  individual  (SN5)  from  HS3:  Example  2. 
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SWEAT  RATES 


Measured  sweat  rate  (SR)  data  from  HS1  and  HS3  were  available.  SR  in  both 
studies  was  calculated  from  difference  in  body  weight  and  associated  weights  (e.g., 
clothing)  before  and  after  work.  Summaries  of  individual  comparisons  between 
measured  and  ICDA  predicted  SR  for  HS1  and  HS3  are  summarized  in  Tables  2  and  3, 
respectively.  For  HS1 ,  the  grand  mean  of  SR  predicted  by  ICDA  (3.23  g/min)  was  not 
statistically  different  from  measured  value  of  SR  (3.13  g/min),  using  a  paired  t-test  (p  < 
0.05). 

Table  2.  Subject  and  group  comparisons  between  measured  Sweat  Rates  (SR)  and 
predicted  Sweat  Rates  (SR)  for  HS1 . 


Subject  ID 

lab  SR  (g/min) 

ICDA  SR  (g/min) 

1 

3.71 

2.39 

2 

3.44 

4.03 

3 

1.57 

2.38 

4 

3.17 

4.35 

5 

3.17 

2.70 

6 

3.97 

3.06 

7 

2.64 

2.90 

8 

2.77 

3.37 

9 

3.71 

3.92 

MEAN 

3.13 

3.23 

STDEV 

0.73 

0.73 

In  comparison  to  HS1,  greater  individual  variability  in  SR  was  observed  under  the  more 
strenuous  heat  conditions  of  HS3.  ICDA  SR  predictions  in  HS3  varied  by  individual: 
similar  SR  predictions  to  measured  SR  (ASR  <  1 .2  g/min)  were  observed  in  three  of 
eight  individuals,  while  A  SR  in  the  rest  of  subjects  varied  between  1 .2  and  6.3  g/min. 

Using  a  paired  t-test  (p 

<  0.05),  mean  values  of  measured  and  ICDA  predicted  SR  were 

not  statistically  different  for  HS3.  However,  despite  the  favorable  results  for  mean 

comparisons, 

individual  variation  for  SR  will  need  to  be  carefully  considered  for  future 

predictions. 

Table  3.  Subject  and  group  comparisons  between  measured  Sweat  Rates  (SR)  and 

predicted  SR  in  HS3 

Subject  ID  lab  SR  (g/min) 

ICDA  SR  (g/min) 

1 

8.83 

3.01 

2 

9.07 

3.64 

3 

7.19 

6.69 

4 

5.28 

10.16 

5 

12.66 

12.46 

6 

11.95 

13.84 

7 

14.33 

8.04 

8 

12.55 

11.41 

MEAN 

10.23 

8.66 

STDEV 

3.12 

4.01 
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DISCUSSION 


The  ICDA  model,  derived  from  a  minimum  number  of  non-invasive  inputs,  was 
designed  primarily  for  screening  physiological  responses  of  deployed  Soldiers  assigned 
for  long-hour  multiple  tasks.  The  initial  analysis  of  this  simple  model  yielded  promising 
results  for  predicting  the  thermal  status  of  Soldiers  during  heat  stress.  The  ICDA 
predictions  for  Tcr  were  within  an  acceptable  range  of  ±1 SD  for  non-heat  acclimated 
volunteers  when  they  were  exposed  to  moderately  stressful  heat  levels  of  27°C/50%RH. 
However,  under  more  stressful  conditions,  the  ICDA  predictions  were  better  for  heat- 
acclimated  individuals  than  non  heat-acclimated  individuals.  Initially  predictions  showed 
Tcr  increasing  faster  than  measured  Tcr  in  heat-acclimated  Soldiers.  However,  by  the 
end  of  exercise,  Tcr  predictions  were  more  accurate.  Under  heat  stress  wearing 
protective  garments,  the  ICDA  predictions  were  close  to  measured  Tcr,  although  error 
estimates  were  larger  when  sample  sizes  were  smaller.  An  option  to  improve  the 
simulation  for  more  diversified  populations  (e.g.,  non-acclimated  persons,  different 
fitness,  and  civilian  occupations)  could  be  added  if  that  is  desirable. 

Predicted  values  for  TSk  were  not  as  accurate  as  Tcr  predictions.  This  may  be 
due  in  part,  to  between  and  within  individual  variability,  and  differences  between  the 
regions  represented  by  different  TSk  measurement  sites.  This  trend  seems  to  be 
especially  true  when  the  environmental  temperature  was  lower  than  measured  TSk-  The 
variability  between  measurement  sites  was  less  when  subjects  wore  protective 
garments  or  when  working  in  a  hot-humid  environment.  Before  more  modeling 
development,  there  should  be  agreement  on  the  best,  most  representative  sites  for 
measurement  of  Tsk  to  accommodate  variability  among  individuals,  clothing, 
environment,  and  operations. 

Two  heat  studies  (HS1 ,  HS3)  allowed  comparison  of  measured  SR  to  the  ICDA 
predicted  SR.  When  the  heat  stress  was  less  strenuous,  the  mean  change  (A)  in  SR 
was  0.1  g/min.  When  the  heat  stress  was  more  strenuous,  the  A  in  SR  was  greater. 
Although  the  mean  predicted  SR  was  not  statistically  different  from  the  measured  SR, 
prediction  errors  increased  when  individual  variation  in  SR  was  increased.  This  may 
indicate  that  the  estimate  of  M,  based  solely  on  HR,  may  not  be  adequate  for  accurate 
physiological  predictions.  For  instance,  predicted  SR  of  two  subjects  in  HS3  was  lower 
than  measured  SR  because  of  the  high  resting  HR.  As  a  result,  a  lower  HR  ratio  was 
used  to  estimate  lower  M  when  Ta  was  constant.  In  these  studies,  resting  HR  was 
determined  at  rest  or  sitting  at  the  beginning  of  the  studies  or  lowest  HR  during  the 
studies.  It  is  possible  that  initial  HRs  were  caused  by  anxiety,  initial  discomfort,  or  other 
physical  or  psychological  factors,  and  were  higher  than  under  “true”  neutral  resting 
conditions.  The  affects  of  using  a  subjects’  “true”  resting  HR  on  the  calculated  value  of 
M  need  to  be  carefully  considered. 

Estimates  of  M  from  HR  obtained  using  non-invasive  methods  are  convenient  for 
Soldiers  who  may  be  required  to  perform  different  tasks  over  many  hours  or  even  days. 
Normally,  Soldiers  deploying  to  regions  known  for  extreme  heat  will  be  acclimatized 
during  training  for  maximizing  their  performance  and  prevention  from  heat-related 
injuries.  As  the  ICDA  was  primarily  developed  for  heat  acclimated  Soldiers,  predictions 
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of  Soldiers’  physiological  status  using  the  current  version  of  ICDA  will  be  more  accurate 
for  acclimated  individuals.  Although  the  present  version  of  ICDA  showed  promising 
results  with  study  data  for  some  physiological  parameters,  further  study,  analysis, 
validation  and  improvements  are  recommended. 
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CONCLUSIONS 


The  human  simulation  computer  program  ICDA  has  been  validated  to  estimate 
Tcr>  TSk,  and  SR  using  HR,  environmental  parameters  and  the  Soldier's  anthropological 
characteristics.  The  predictions  were  compared  to  measured  values  from  three 
laboratory  heat  studies.  Overall  the  simulation  predictions  of  Tcr  agreed  well  with  the 
measured  values  for  heat  acclimatized  Soldiers.  The  agreement  for  TSk  and  SR  was 
less  accurate  than  Tcr  and  varied  with  clothing,  environment,  activity  and  individual 
differences,  and  should  be  evaluated  further.  Additional  future  experience  with  ICDA 
during  field  training  situations  will  lead  to  further  assessment,  confidence,  and  possible 
improvements  in  application  and  hardware. 

The  ICDA  real  time  physiological  monitoring  system  was  created  for  the  narrow 
specialized  application  of  trained  heat  acclimated  Soldiers  ready  for  deployment  to  a  hot 
environment.  However  the  system  can  be  adapted  for  more  variation  in  fitness,  heat 
acclimation,  body  mass,  age  and  other  characteristics  for  application  to  more  diverse 
groups  of  population  and  environment. 


-22  - 


REFERENCES 


Amos,  D.,  R.  Hansen,  W.  Lau,  J.T.  Michalski.  Physiological  and  cognitive  performance 
of  Soldiers  conducting  routine  patrol  and  reconnaissance  operations  in  the  tropics.  Mil. 
Med.  165:  961-966,  2000. 

Berglund,  L.G.  Heart  rate  as  an  indicator  of  metabolic  rate  in  hot  environments.  30th 
Annual  Conference  on  Engineering  in  Medicine  and  Biology,  19:  274,  1977. 

Bilski,  T.R.,  B.C.  Baker,  J.R.  Grove,  R.  P.  Hinks.  Battlefield  casualties  treated  at  Camp 
Rhino,  Afghanistan:  Lessons  learned.  J.  Trauma.  54:814-822,2003. 

Fanger,  P.O.  Thermal  Comfort.  New  York:  MacGraw  Hill.  1972. 

Haslam,  R.  A.,  and  K.  C.  Parsons.  An  evaluation  of  computer-based  models  that  predict 
human  responses  to  the  thermal  environment.  ASHRAE  Trans.  94:  1342-1360,  1988. 

Gagge,  A.  P.,  A.  P.  Fobelets,  L.G.  Berglund.  A  standard  predictive  index  of  human 
response  to  the  thermal  environment.  ASHRAE  Trans.  92(1),  1986. 

Gawande,  A.  Casualties  of  war  -  Military  care  for  the  wounded  from  Iraq  and 
Afghanistan.  N.Engi.J.Med.  351:2471-2475,2004. 

Kraning,  K.  K.,  R.  R.  Gonzalez.  Physiological  consequences  of  intermittent  exercise 
during  compensable  and  uncompensable  heat  stress.  J.  Appl.  Physiol.  71:2138-2145, 
1991. 

Kraning,  K.  K.,  R.  R.  Gonzalez.  A  mechanistic  computer  simulation  of  human  work  in 
heat  that  accounts  for  physical  and  physiological  effects  of  clothing,  aerobic  fitness,  and 
progressive  dehydration.  J.  Therm.  Biol.  22,331-342,  1997. 

Latzka,  W.  A.,  M.  N.  Sawka,  S.  J.  Montain,  G.  S.  Skrinar,  R.  A.  Fielding,  R.  P.  Matott,  K. 
B.  Pandolf.  Hyperhydration:  tolerance  and  cardiovascular  effects  during 
uncompensable  exercise-heat  stress.  J.  Appl.  Physiol.  84,1858-1864,  1998. 

Montain,  S.  J.,  S.  N.  Cheuvront,  L.  G.  Berglund,  A.  J  Young,  M.  N.  Sawka,  R.  Carter. 
Prediction  and  validation  of  Soldier  water  requirements.  Natick,  MA:  USARIEM. 

Human  Research  Protocol  H03-14,  2005. 

Moran  D.  S.,  A.  Shitzer,  K.  B.  Pandolf.  A  physiological  strain  index  to  evaluate  heat  stress. 
Am.  J.  Physiol.  275:R129-R134,  1988. 

Pandolf,  KB,  and  Goldman,  RF.  Convergence  of  skin  and  rectal  temperatures  as  a 
criterion  for  heat  tolerance.  Aviat.  Space  Environ.  Med.  49:1095-1 101,  1978. 


-23- 


Ramanathan,  N.  L.  A  new  weighting  system  for  mean  surface  temperature  of  the  human 
body.  J.  Appl.  Physiol.  19:  531-533,  1964. 

Santee,  W.  R.,  L.  G.  Berglund,  A.  Cardello,  C.  A.  Winterhalter,  T.  L.  Endursick. 
Physiolgoical  assessment  of  volunteers  wearing  Battle  Dress  Uniforms  (BDU)  during 
intermittent  exercise.  Natick,  MA:  USARIEM.  Technical  Report  T05,  2005. 


-24- 


APPENDIX  A 


A  documented  listing  of  the  basic  ICDA  model  described  in  the  introduction  and 
used  in  the  verification  follows.  It  is  slightly  different  from  the  real  time  simulator  on  the 
Soldier  in  that  measured  input  is  entered  through  a  keyboard  rather  than  directly  from 
the  sensor  system.  Also  the  specified  measurement  interval  is  not  driven  by  a  clock  in 
this  case  so  the  measured  values  can  be  keyboard  entered  at  operators  pace. 

Following  the  program  instructions  an  example  of  the  output  is  given. 

/*  Origin:  BBMD  @  USARIEM 
Authors:  Berglund  and  Yokota 
Date:  6/30/05 

File:  IC-DA-C-RT_Beta__steppedNonClockedlnputs6_30_05 

This  simulates  the  real  time  thermal  physiological  responses(Tc,Tsk,waterLoss,skin  wettedess,  PSIc) 
from  real  time  inputs  of  measured  HR,  environmental  parameters(Ta,MRT,RH,V)  with  constant  personal 
properties(ht,wt)  and  clothing. 

The  interval  (in  seconds)  between  measurements  is  adjustable.  Enter  0  for  the  default  interval  of  60 
seconds.  In  this  program  version  the  input  interval  is  not  controlled  by  a  clock. 

Measurements  values  are  entered  through  key  board  at  these  intervals.  The  measurement  interval  is  also 
the  model's  integration  interval. 

Metabolism  is  expressed  by  dimenionless  relative  metabolism  term  (met)  where  met=actual 
metabolism/resting  metabolism.  Met  is  estimated  in  this  program  from  heart  rate(HR)  and  air 
temperature(Ta).  intrinsic  clo  is  thermal  resistance  from  skin  to  outer  layer  of  clothing  in  clo  units. 

MRT=  mean  radiant  temperature.  V=air  speed.  7 

#include  <assert.h> 

#include  <stdio.h> 

#include  <math.h> 

#include  <ctype.h> 

#define  TTSK  33.7  //°C 
#define  TTCR  36.8  //°C 

#define  TTBM  36.49  //  TTBM=.9*TTCR+.1*TTSK  =36.49 
#define  CSW  170.  //  g/(h  mA2  °C) 

#define  CDIL  50.  //  CDIL  200  super  athlete,  50  average  person 

#define  CSTR  .5  //  1/°C 

#define  SKBFN  6.3  //  liters/(h  mA2) 

#define  Skbfmax  90.  //  conservative  could  be  higher  for  fit  person 
#define  Skbfmin  2.  //Liters/(h  mA2) 

#define  CMIN  5.28  //  w/(mA2  °C) 

#define  CB  1.163  //wh/(L°C) 

#define  WpMet  58.1  //  w/(m2*met) 

#define  LR  2.2  //Lewis  relation 

void  OutputHeader(); 
double  SatVapPres(doub!e  T); 

double  Clot(int  code,  double  V);  //total  clo  from  skin  to  environment 
double  VRes(int  code,  double  V);  //  total  vapor  resistance  from  skin  to  environment 

double  Convection(double  met,  double  V);  //convective  heat  transfer  coefficient 
double  Respiration(double  Rm,  double  Ta,  double  Pa, double*  Wres);  //respiratory  heat  loss 
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double  SkinBloodFlow(double  Tc,  double  Tsk); 

double  DuBoisSkinArea(double  Wt,  double  Ht);  //surface  area  of  body 

double  MetT(double  HR, double  HRrest, double  T); 

double  Alpha(double  Skbf);  //  fraction  of  body  that  is  skin 

double  Sweat(double  Tc,  double  Tsk, double  alpha); 

double  Core(double  Tc, double  Tsk, double  met, double  Res, double*  hfcsk,  double  skbf, double*  Rm); 

double  Skin(double  wet, double  Dry, double  Emax, double*  hfcsk); 

double  MoranPSI(double  Tc, double  Tco, double  HR, double  HRrest);  //heat  stress  index 

int  main  (void){ 

int  step,  HrTimeStep,code;double  Ta,Tanew,MRT,MRTnew,RH,RHnew,Pa,Psk,Tcl,To,V,Vkm,\/new; 

double  hc,he,hr,met,Ht,Wt,Adu,Wres,  drip,DT,DTs,ExTime; 

double  wet,Tsk,Tsko,Tc,Tco,Rm,cio,FACL,Dry,Fpcl,PSIc; 

double  Emax,Esk, Edit, Res, Skbf, HRrest, HR, HRnew, alpha,  Regsw.Ersw; 

double  heatFlowCoreToSkin,HSCR,HSSK,TCCR,TCSK,Rcl,Rbound,Rclt,Rpcl,  Rpbound,Rpclt; 

double  time, TIM,  Duration, DurationM.CumulativeWaterLossPsqm; 

doubleCumulativeWaterLoss,CWLL,RateOfTotalWaterLossPsqm,  RateOfTotalWaterLoss,Wdif; 

TIM=0;  step=0;  time  =0;  HR=0;  HRrest-0;  DTs=0;ExTime=0;  HrTimeStep=0; 

printf("Enter  interval  between  measurements  in  seconds(0  for  default  of  60  seconds):DTs  "); 

scanf("%lf",&DTs); 

if  (DTs<=0) 

DTs=60; 

DT=DTs/(60*60);  //DT  in  hrs 

printf("Enter  intrinsic  clo  or  -1  if  unknown:  clo  "); 

scanf("%lf',&clo); 

code=0; 

if  (clo<0){ 

printf("  Codes  of  manikin  tested  clothing  types  available  for  simulation:  \n  1  HWBDU  \n  2 
APFU(SHTS+T)  \n  3  Protective  Garment(MOPP4)  \n  Enter  Code:  code  "); 
scant  ("%d",&code);  } 

printffEnter  subject's  weight(kg),height(m)(or  0's  if  unknown):  Wt  Ht  ");/*  inserts  default  values  of 
70kg,  1.72m  7 
scanf("%lf%lf',&Wt,&Ht); 
if  (W  t<— 0) 

Wt=72.8; 
if  (Ht<=0) 

Ht=1 .72; 

printf("Enter  subject's  resting  HR)(or  0's  if  unknown):HRrest  ");/*  inserts  default  values  of  70  if  HRrest 
are  zero  7 
scanf("%lf",&HRrest); 
if  (HRrest<=0) 

HRrest=70; 

Adu=DuBoisSkinArea(Wt,Ht);  //mA2 

printf("Enter  initial  physiology  (or  0's  if  unknown):  Tsk  Tc  ");/*  inserts  default  values  if  Tsk, Ter  are 
zero's  7 

scanf(''%lf%lf",&Tsk,&Tc); 
if  (Tsk  <=0){ 

Tsk=33;  Tsko=Tsk; } 

else 

Tsko=Tsk; 
if  (Tc  <=0){ 

Tc=36.9;  Tco=Tc;  } 

else 

Tco=Tc; 
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wet=.06; 

printf("Enter  duration(min):  ");  //may  wish  to  delete  duration  for  real  applications 

scanf("%lf",&DurationM); 

Duration=DurationM/60.;  //hr 

printffEnter  HR  Ta(C)  MRT(C)  RH(percent)  V(km/h) :  "); 

scanf("%lf%lf%lf%lf%lf",&HR,&Ta,&MRT,&RH,&Vkm); 

V  =Vkm*1 000/(60*60);  //  convert  km/h  to  m/s 
met=MetT(HR,HRrest,Ta);  //from  Berglund  30th  ACEMB.1977 
OutputHeader(); 

Pa=RH*SatVapPres(Ta)/100;  //  vapor  press  Torr 
Skbf=SkinBloodFlow(Tc,Tsk);  //Liters/(h  mA2) 
alpha=Alpha(Skbf);  //skin  fraction  of  body  =Wtsk/Wt 
Regsw=Sweat(Tc,Tsk, alpha);  //  g/(h  mA2) 

CumulativeWaterLossPsqm=0; 

CumulativeWaterLoss=0; 
drip=0;  CWLL=0; 

he  =  Convection(met,V);  //convection 
hr— 4.5;  //radiation  watts/(mA2  K) 

To=(hc*Ta  +  hr*MRT)/(hc+hr); 
if  (code==0)  { 

FACL=1 ,0+.2*c!o;  //  surface  area  of  clothing  relative  to  Adu 
Rcl=.155*clo;  //  clothing  thermal  resistance  mA2  C/watts 
Rbound=1/((hr+hc)*FACL);  //  boundary  layer  thermal  resistance  mA2  C/watts 
Rclt=Rcl+Rbound; 

Rpck  ,153181*clo;  } 

else  { 

FACL=1 .2;  //guess  with  clo=1 

Rclt=.155*Clot(code,V); 

Rbound=1/((hr+hc)*FACL);  //  boundary  layer  thermal  resistance  mA2  C/watts 
Rcl=Rclt-Rbound; 

clo=Rcl/.155;  //estimate  of  uniform's  intrinsic  clo 
FACL=1+.2*clo; } 

Dry=(Tsk-To)/(Rclt); 

Tel  =  To  +Dry/(FACL*(hr+hc)); 

PSic=MoranPSI(  Tc,Tco,  HR,  HRrest); 

printf("%5.1f  %5.1f  %5.1f%4.0f  %5.2f  %5.2f  %5.2f  %5.2f  %5.2f  %5.2f%4.0f  %5.3f  %5.2f  \n 

",ExTime,Ta,MRT,RH,Vkm,Tc,Tsk,met,Regsw/60*Adu,wet,HR,CWLL,PSIc); 

time+=Duration; 

while  (TiM<=time)  {  //  thermo-physioiogy  ioop 
//  dry  and  evaporative  heat  transfer 

Skbf=SkinBloodFlow(Tc,Tsk);  //Liters/(h  mA2) 

Psk=SatVapPres(Tsk);  //  vapor  press  Torr  (mmHg) 

Rm=  WpMet*met;  //  watts/(mA2) 

he  =  Convection(met,V);  //  convection  */ 

he=2.2*hc;  //evaporation  watts/(mA2  Torr) 

hr=4*.725*(5.67E-08)*pow(((Tcl+To)/2+273),3);  //corrected  hr 

To=(hc*Ta  +  hr*MRT)/(hc+hr); 

if  (code==0)  { 

Fpcl=1/(1  +.1 531 81  *he*FACL*clo);  //Berglund  1981,  for  IL=.45 
Rpbound=1  /(he*FACL); 

Rpclt=Rpcl+Rpbound; 

Rbound=1/((hc+hr)*FACL); 

Rclt=Rbound+Rcl;  } 
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else  { 

Rpclt=VRes(code,V); } 
Emax=(Psk-Pa)/Rpclt; 

Dry=(T  sk-T  o)/(Rclt); 

Tel  =  To  +Dry/(FACL*(hr+hc)); 


//  Thermal  Physiology 

Res  =  Respiration(Rm,Ta,Pa,&Wres); 
alpha=Alpha(Skbf); 

Regsw=Sweat(T  c,Tsk, alpha); 
Ersw=.68*Regsw; 

Edit  =  (1-wet)*.06*Emax; 
Wdif=Edif/.68; 


Esk=(Wdif+Regsw)*.68; 
wet=Esk/Emax; 
if  (wet>=1 ){ 
wet=1 ; 


//  g/(h  mA2) 

//  watts/mA2 

//  diffusion  through  dry  skin 
//  g/(h  mA2) 


//  energy  lost  by  evaporation  from  and  diffusion  through  skin. 


drip=(Wdif+Regsw)-(Emax/.68); }  //  wasted  sweat  driping  from  skin  g/(hmA2) 

else  { 
drip— 0; } 

Esk=(.06+.94*wet)*Emax; 

RateOfTotalWaterLossPsqm=Regsw+Wres+Wdif;  //  g/(h  mA2) 

RateQfT otalWaterLoss=RateOfT otalWaterLossPsqm*Adu;  //g/h  mA2) 
CumulativeWaterLossPsqm=CumulativeWaterLossPsqm+RateOfTotalWaterLossPsqm*DT; 
CumulativeWaterLoss=CumulativeWaterLoss+RateOfTotalWaterLoss*DT; 
CWLL=CumulativeWaterLoss/1 000;  //Liters 

PSIc=MoranPSI(  Tc.Tco,  HR,  HRrest); 

HSCR  =  Core(Tc,Tsk,met,Res,&heatFlowCoreToSkin,Skbf,&Rm); 

HSSK  =  Skin(wet,Dry,Emax,&heatFlowCoreToSkin); 


/*  thermal  capacity  V 

TCCR=.97*(1  -alpha)*Wt; 
TCSK=.97*alpha*Wt; 


/*  stepwise  integration  V 

Tc  +=HSCR*Adu/TCCR*DT; 

Tsk  +=HSSK*Adu/TCSK*DT; 

TIM  +=DT; 

ExTime+=DT*60;  //min 

printf("Enter  HR  Ta  MRT  RH  V  or  0  for  any  unchanged  parameter:  "); 
scanf("%lf%lf%lf%lf%lf",&HRnew,&Tanew,&MRTnew,&RHnew,&Vnew); 
if  (HRnew>0) 

HR=HRnew; 

If  (Tanew>0) 

Ta=Tanew; 

If  (MRTnew>0) 

MRT=MRTnew; 

If  (RHnew>0) 

RH=RHnew; 

If  (Vnew>0)  { 

Vkm=Vnew;  //km/h 

V=Vkm*1 000/(60*60); }  //m/s 

met=MetT(HR,  HRrest, Ta); 

printf(“%5.0f  %5.1f  %5.1f%4.0f  %5.2f  %5.2f  %5.2f  %5.2f  %5.2f  %5.2f%4.0f%5.3f%5.2f\n 
",ExTime,Ta,MRT,RH,Vkm,Tc,Tsk,met,Regsw/60*Adu,wet,HR,CWLL,PSIc); 
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return  0  ; 


void  OutputHeader()  { 

printf("\nmin  Ta  Tr  RH  V  Tc  Tsk  met  sw  wet  HR  CWLoss  PSIc\n"); 
printf("  C  C  km/h  C  C  g/min  bpm  L\n"); 

} 

double  SatVapPres(double  T)  { 
double  Pst; 

Pst  =  exp(18.6686-(4030. 1 83/(T+235.)));  //mmhg 

return  Pst; 


double  Clot(int  code,  double  V)  { 
double  A1=0,B=0; 
switch  (code)  { 

case  1 :  //HWBDU 

A1=1.08;  B=-0.27;  break; 
case  2:  //APFU(SHTS+T) 

A1=0.51;  B=-0.39;  break; 
case  3:  //Protective  Garment(MOPP4) 

A1=l  .7;  B=-0.16;  break; 

default:printf("lllegal  clothing  code, must  be  an  integer  between  1  and  3  ");break; } 
return  A1*pow(V,B); 

} 

double  VRes(int  code,  double  V)  { 
double  Ai1=0,Bi=0; 
switch  (code)  { 

easel:  //HWBDU 

Ail  =0.47;  Bi=0.41;  break; 
case  2:  //APFU(SHTS+T) 

Ail  =1.41;  Bi=0.58;  break; 
case  3:  //Protective  Garment(MOPP4) 

Ail  =0.1 7;  Bi=0.31;  break; 

default:printf("lllegal  clothing  code, must  be  an  integer  between  1  and  3  ");  break; } 
return  0.155/(LR*  Ai1*pow(V,Bi)); 


double  Convection(double  met, double  V)  ( 

double  CHCA,  CHCV,  CHCmin,  he;  CIHCmin  =  3.0; 
if  (met>=1 .1)  { 

CHCA  =  5.66*pow((met  -  0.85), 0.39);  //  he  due  to  activity 

else  { 

CHCA  =  5.66*pow((1 .  -  0.85), 0.39);  }  //  he  due  to  activity 

} 

CHCV  =  8.6*pow(V,0.53);  //  he  due  to  air  speed  V  in  m/s  7 

if  (CHCV  >=  CHCmin)  ; 

else 

(CHCV  =  CHCmin); 

if  (CHCV  >=  CHCA) 
he  =  CHCV; 

else 

he  =  CHCA; 
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return  he; 


} 

double  Respiration(double  Rm,  double  Ta,  double  Pa,double*Wres){ 
double  Res,Eres,Cres; 

Eres  =  0.0023*Rm*(44.-Pa);  //  watts/m2 

Cres  =  0.001 4*Rm*(34.-Ta);  //  watts/m2 

*Wres=Eres/.68;  //  g/(h  mA2)  7 

Res  =  Eres  +  Cres; 
return  Res; 

} 

double  DuBoisSkinArea(double  Wt,  double  Ht)  { 
double  Adu;  /*  mA2  7 

Adu=0.202*pow(Wt,0.425)*pow(Ht, 0.725);  //  Wt  weight  in  kg,  Ht  height  in  m. 
return  Adu; 

} 


double  MetT(double  HR, double  HRrest, double  T)  {  //  Berglund,30th  ACEMB  1977 

double  met,HRratio;  /*  applicable:20<=T<=40,  0.6clo,  v~=1.25m/s,  1.2<+HRratio<=2.1,Tdp<=20C7 
HRratio^HR/HRrest; 
if  (HRratio>=1.2)  { 
if  (T>=20) 

met=0.68+4.69*(HRratio  -1)-  0.052*(HRratio-1)*(T-20); 

else 

met=0.68+4.69*(HRratio  -1);  } 

else  { 

met=HRratio;  } 
return  met; 


double  Alpha(double  Skbf)  { 
double  alpha; 

alpha=0.04177+.74518/(Skbf+0.585417); 
return  alpha; 


double  Sweat(doubie  Tc,  double  Tsk,  double  alpha)  { 
double  regsw,  Tmb; 
regsw=0; 

Tmb=(1-alpha)*Tc  +  alpha‘Tsk; 
if  ((Tmb>TTBM)&&(Tsk>TTSK)) 

regsw=CSW*(Tmb-TTBM)*exp((Tsk-TTSK)/10.7); 
else  if  ((Tmb>TTBM)&&(Tsk<=TTSK)) 
regsw=CSW*(T  mb-TTBM); 

if  (regsw>667)  //  max  sweat  rate  limit 

regsw=667;  //  regsw_max=667g/(h  mA2)=1 1 ,1g/(min  mA2)!=20g/(min  mA2) 

return  regsw;  //g/(h  mA2) 


double  SkinBloodFlow(double  Tc,  double  Tsk)  { 
double  Colds=0; 
double  Skbf,  WarmC=0; 
if  (Tsk<TTSK) 

Colds=TTSK-Tsk; 
if  (ToTTCR) 
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Warm  C=T  c-TT  C  R ; 

Skbf=(SKBFN+CDIL*WarmC)/(1  +CSTR*Colds);  //  Liters/(h  mA2) 

if  (Skbf>Skbfmax) 

Skbf=  Skbfmax; 
if  (Skbf<Skbfmin) 

Skbf=  Skbfmin; 

return  Skbf;  //  L/(h  mA2) 


double  Core(double  Tc, double  Tsk, double  met, double  Res,double*heatFlowCoreToSkin,  double  Skbf, 
double*Rm)  { 

double  Rmet,Hfcrsk,FISCR; 

Hfcrsk=(CMIN+CB*Skbf)*(Tc-Tsk);  //  watts/mA2 

Rmet  =  58.2*met ;  //  metabolic  heat  produced  watts/mA2. 

HSCR=Rmet-Hfcrsk-Res;  //  rate  of  heat  storage  in  core  watts/mA2 

‘heatFlowCoreT  oSkin=Hfcrsk; 

*Rm=Rmet; 
return  FISCR; 


double  Skin(double  wet, double  Dry, double  Emax,double*heatFlowCoreToSkin  )  { 
double  Esw,Ediff,Esk,FISSK; 

Esw=wet*Emax; 

Ediff=.06*(1-wet)*Emax; 

Esk=Esw+Ediff; 

HSSK=*heatFlowCoreToSkin-Dry-Esk;  //  rate  of  heat  storage  in  skin  watts/mA2 
return  FISSK; 

} 

double  MoranPSI(double  Tc, double  Tco, double  FIR,  double  HRrest){ 
double  PSI; 

PSI=5*((Tc-Tco)/(39.5-Tco)+(HR-HRrest)/(180-HRrest)); 
return  PSI; 

} 

/*  This  is  the  End  7 


/*  output  is  pasted  below 

Enter  interval  between  measurements  in  seconds(0  for  default  of  60  seconds):DTs  0 

Enter  intrinsic  clo  or  -1  if  unknown:  clo  -1 

Codes  of  manikin  tested  clothing  types  available  for  simulation: 

1  FiWBDU 

2  APFU(SHTS+T) 

3  Protective  Garment(MOPP4) 

Enter  Code:  code  2 

Enter  subject's  weight(kg),height(m)(or  0's  if  unknown):  Wt  Ht  0  0 
Enter  subject's  resting  FIR)(or  0's  if  unknown):FIRrest  0 
Enter  initial  physiology  (or  0's  if  unknown):  Tsk  Tc  0  0 
Enter  duration(min):  10 

Enter  FIR  Ta(C)  MRT(C)  RFI(percent)  V(km/h)  :  75  27  27  75  4 

min  Ta  Tr  RFI  V  Tc  Tsk  met  sw  wet  FIR  CWLoss  PSIc 
C  C  km/h  C  C  g/min  bpm  L 

0.0  27.0  27.0  75  4.00  36.90  33.00  1.07  0.00  0.06  75  0.000  0.23 
Enter  HR  Ta  MRT  RH  V  or  0  for  any  unchanged  parameter:  90  0  0  0  0 
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1  27.0  27.0  75  4.00  36.90  32.78  1.92  0.00  0.06  90  0.001  0.23 
Enter  HR  Ta  MRT  RH  V  or  0  for  any  unchanged  parameter:  109  0  0  0  0 

2  27.0  27.0  75  4.00  36.92  32.58  3.09  0.00  0.06  109  0.002  0.91 
Enter  HR  Ta  MRT  RH  V  or  0  for  any  unchanged  parameter:  107  0  0  0  0 

3  27.0  27.0  75  4.00  36.98  32.41  2.97  0.00  0.06  107  0.004  1.82 
Enter  HR  Ta  MRT  RH  V  or  0  for  any  unchanged  parameter:  108  0  0  0  0 

4  27.0  27.0  75  4.00  37.02  32.27  3.03  0.00  0.06  108  0.005  1.83 
Enter  HR  Ta  MRT  RH  V  or  0  for  any  unchanged  parameter:  0  0  0  0  0 

5  27.0  27.0  75  4.00  37.06  32.16  3.03  0.00  0.06  108  0.006  1.96 
Enter  HR  Ta  MRT  RH  V  or  0  for  any  unchanged  parameter:  0  0  0  0  0 

6  27.0  27.0  75  4.00  37.10  32.06  3.03  0.26  0.07  108  0.008  2.04 
Enter  HR  Ta  MRT  RH  V  or  0  for  any  unchanged  parameter:  0  20  20  0  0 

7  20.0  20.0  75  4.00  37.13  31.98  3.23  0.48  0.09  108  0.010  2.11 
Enter  HR  Ta  MRT  RH  V  or  0  for  any  unchanged  parameter:  0  0  0  0  0 

8  20.0  20.0  75  4.00  37.17  31.51  3.23  0.68  0.10  108  0.012  2.17 
Enter  HR  Ta  MRT  RH  V  or  0  for  any  unchanged  parameter:  0  0  0  0  0 

9  20.0  20.0  75  4.00  37.20  31.12  3.23  0.52  0.09  108  0.013  2.24 
Enter  HR  Ta  MRT  RH  V  or  0  for  any  unchanged  parameter:  0  0  0  0  0 

10  20.0  20.0  75  4.00  37.23  30.81  3.23  0.43  0.09  108  0.015  2.30 
Press  any  key  to  continue 

7 
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APPENDIX  B 


The  mean  HR  and  ±1SD  from  each  study  are  displayed, 
a.  Heat  Study  1  (27°C,  75%;  HWBDU;  n=9) 
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b.  Heat  Study  2  (dayl  -  non  acclimation  state;  49°C,  18%;  APFU;  n<7) 
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APPENDIX  B  cont. 


c.  Heat  Study  2  (day2  -  acclimation  state;  49°C,  18%;  APFU;  n=5) 
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d.  Heat  Study  3  (35°C,  45%;  protective  garment;  n<8) 
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